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ABSTRACT 

Motivated by the possibility that a coronagraph will be put on WFIRST/AFTA, we explore the 
direct detectability of extrasolar giant planets (EGPs) in the optical. We quantify a planet’s detectabil¬ 
ity by the fraction of its orbit for which it is in an observable configuration (/obs)- Using a suite of 
Monte Carlo experiments, we study the dependence of /obs upon the inner working angle (IWA) and 
minimum achievable contrast (Cmin) of the direct-imaging observatory; the planet’s phase function, 
geometric albedo, single-scattering albedo, radius, and distance from Earth; and the semi-major axis 
distribution of EGPs. We calculate phase functions for a given geometric or single-scattering albedo, 
assuming various scattering mechanisms. We find that the Lambertian phase function can predict 
significantly larger /obs’s with respect to the more realistic Rayleigh phase function. Eor observations 
made with WEIRST/AETA’s baseline capabilities (Cmin ^ 10“^, IWA ^ 0.2"), Jupiter-like planets 
orbiting stars within 10, 30, and 50 parsecs of Earth have volume-averaged observability fractions 
of ^12%, 3%, and 0.5%, respectively. At 10 parsecs, such observations yield /obs > 1% foi* low- to 
modest-eccentricity planets with semi-major axes in the range ^2 — 10 AU. If Cmin = 10“^^, this 
range extends to ^35 AU. We find that, in all but the most optimistic configurations, the probability 
for detection in a blind search is low (< 5%). However, with orbital parameter constraints from long¬ 
term radial-velocity campaigns and Gaia astrometry, the tools we develop in this work can be used 
to determine both the most promising systems to target and when to observe them. 

Subject headings: planetary systems - planets and satellites: general 


1. INTRODUCTION 

In the burgeoning field of exoplanetary science, the sci¬ 
entific return from direct imaging has so far been meager 
when compared to the very successful radial-velocity and 
transit methods. This delay in progress has been pri¬ 
marily due to the technological challenge of performing 
high-contrast observations at small angular separations 
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enheimer & Hinkley 2009 Traub & Oppenheimer 


kbr a Jupiter twin, the plant/star hux ratio is ex- 
pected to be ^10“^ in the mid-infrared and r\j 10 ^ in the 
optical, and if the planet is in close proximity to Earth 
(< 30 pc), these contrasts must be achieved at angular 
separations of a few tenths of arcseconds. Making the 
problem harder still, the direct detectability of a given 
planet will vary with time, depending upon its Keplerian 
orbital elements and atmospheric properties — which will 
depend upo n the planet’s age, mass, orbital distance, and 
parent star ([Sudarsky et aL]|2005 _ Kane fc Gelino 2010 
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TJespite these challenges, recent progress has been 
made via high-contrast observations in the near-infrared 
of nearby, young extrasolar giant planets (EGPs)/brown 
dwarfs at very wide orbital separations (^10 — 200 AU). 
Such substellar companions are self-luminous and pro¬ 
vide contrasts and angular separations that are far more 
favorable than the case of a nearby Jupiter twin. In¬ 
deed, there have been a number of notable discoveries 
and partial character izations, inc luding 2MASS 1207b 


(Ghauvin et al.||2004), AB Pic b ([^ 
GCj Lup b ( [P^euhauser et al. 2ofef 7 
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ge et al.||2009E HR 8799bcde (iMa rois et al.||2008[ 
|Barman et al.|201l|) , and G J 504b (|Kuzuhara et al 
2013]). In the coming years, a plethora of high-contrast 


instrument s such as t he ground-based G PI ([Ma cintosh 
et al.|200"8|), SP HERE (jBeuzit et al.|2008), LBTl (Iskrut- 
skie e t al.IzOloj), ScEx AU / Charis / Hid AO (ISuzuki et al. 
2010|) ancl the space-based NIRCam and MIRT 


ScEx Au / Charis / Hid AO 
and the space- based NIRCam and MIRI on J WST 
( jPeming et al.|20(i9 ) will almost certainly produce a rich 
sample of direct lymiag ed and characterized young EGPs 
(Beichman et alj[2010 ). 


All of the above high-contrast imaging efforts are fo¬ 
cused on measurements in the infrared. However, opti¬ 
cal measurements of albedo spectra, phase curves, and 
polarization—which have been used to infer chemical 
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the solar system ([Hansen & Hovenier 
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at wide orbital distances, a population of planets that 
has been inaccessible to infrared measurements. Direct 
imaging and characterization of cool, mature EGPs will 
require a space-based high-contrast instrument such as 
the anticipated optical imaging and spectrophotometric 
coronagraph that is being designed for the Wide-Eield 
Infrared Survey Telescope (WEIRST)/Astrophysics Eo- 
cused Telescope Assets (AETA) mission. If realized, this 
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instrument will achieve planet/star contrast ratios better 
than ^10“^ and an inner worki ng angle near or better 
than ^0.2" (Spergel et ah 2015|), well within the neces¬ 
sary capability to detect a Jupiter twin located ^20 pc 
form EartlQ Therefore, studies of the observational sig¬ 
natures and direct detectability of reflected light from 
EGPs are particularly timely and will aid in the design 
and optimization of WEIRST/AETA or a similar space- 
based effort. 

In this work, we present an exploration of the di¬ 
rect detectability of cool EGPs at visible wavelengths 
(^0.4—1.0 /im). We quantify a planet’s direct detectabil¬ 
ity by the fraction of its orbit for which it is in an observ¬ 
able configuration (its observability fraction /obs)- Using 
a suite of Monte Carlo (MG) experiments, we explore 
the dependence of /obs upon various input parameters 
such as the inner working angle (IWA), minimum achiev¬ 
able contrast (Cmin), planetary albedo and radius, dis¬ 
tance from Earth, and the semi-major axis distribution 
of EGPs. We assume cloud-free, homogeneous atmo¬ 
spheres and calculate phase functions for a given geomet¬ 
ric or single-scattering albedo, assuming different scatter¬ 
ing mechanisms. Although partly motivated by the pos¬ 
sibility of an optical coronagraph on WEIRST/AETA, 
our study is relevant to any effort to directly detect ex¬ 
oplanets in the optical. 

In ^we describe how we model reflected-light observa¬ 
tions of EGPs. We tabulate and provide analytic fits to 
the geometric albedo as a function of single-scattering 
albedo for scalar and vector Rayleigh, Lambert, and 
isotropic scattering. In addition, we present represen¬ 
tative albedo spectra for homogeneous atmospheres of 
cool EGPs. This section is supported by an appendix, 
which reviews the equations that connect the Keplerian 
orbital elements to the observed phase angle and time. 
In ^we present illustrative calculations of the variation 
of planet/star flux ratios and angular separations with 
orbital parameters, which ultimately determines an exo¬ 
planet’s direct detectability. In ^we describe our MG 
experiments and discuss the assumptions we make about 
the orbital parameter distributio ns of EGPs. 

Our results are given in ^ In §5.1|we present a para¬ 


metric study of the mean observability fraction, (/obs), 
for many MG experiment s, wh ich are run under different 
sets of assumptions. In §5.2| we highlight the fact that 
(/obs) for a population of planets is proportional to the 
probability that a blind search for such planets will yield 
a detection. Equivalently, (/obs) can be interpreted as 
the probability f or d etection, assuming all stars have one 
such planet. In §5.3| we briefly consider the optical signa¬ 
ture of thermal emission from young EGPs. The paper 
is brought to close in ^with a summary of our results. 

Throughout this paper, we define a “Jupiter-like” 
planet to be an EGP with the radius of Jupiter, a geomet¬ 
ric albedo of 0.5, and a Rayleigh scattering atmosphere. 
Although these atmospheric properties are not consis¬ 
tent with Jupiter at all wavelengths, near ^0.565 /im, 
which is at the center of one of the WEIRST/AETA 
bandpasses , they provide a d ecent approximation to ob¬ 
servations ( Karkoschka||1994 ). 


^ For a summary of the diagnostic potential for cold exoplanet 
characteri zation using th e anticipated WFIRS T/AFTA corona- 
graph, see [Burrows] ( [20T4| , [Marley et ah] ( |2014| ) , and |Hu| ( |2014| ). 


2. MODELING REFLECTED-LIGHT OBSERVATIONS OF 
GIANT EXOPLANETS 

To directly detect an extrasolar planet, its dim light 
must be separated from the overwhelming glare of its 
bright parent star. Eor a given high-contrast imaging ob¬ 
servation, this requires that two conditions be simultane¬ 
ously satisfied: 1) the planet/star contrast ratio {Fp/F^) 
must be above the minimum achievable contrast (Cmin), 
and 2) the planet/star angular separation (J0) must be 
both larger than the inner working angle (IWA) and 
smaller than the outer working angle (OWA). As a fur¬ 
ther complication, Cmin will typically be a function of 
and both of these parameters will depend upon the wave¬ 
length of the observation. Thus, the direct detectability 
of an exoplanet ultimately depends upon its location in 
{Fp/F^) — 5(j) parameter space, which is a complicated 
function of time, orbital and atmospheric parameters, 
and the high-contrast imaging technology used to make 
the observation. In this section, we describe the methods 
we employ to predict the direct detectability of EGPs as 
they move through (Fp/F^) — 5(j) parameter space. 


2.1. Planet/Star Contrast Ratios: Albedos and Phase 
Funetions 

At optical wavelengths (^0.4 —1.0 fim)^ the flux from a 
giant exoplanet is dominated by reflected starlight, pro¬ 
vided the planet is sufficiently old and of low eno ugh mass 
that its thermal emission is negligible (see ^5.3 for a dis¬ 
cussion of thermal emission from young EGPs). Obser- 
vationally, the wavelength- and phase-dependence of this 
reflection is encompassed by the geometric albedo spec¬ 
trum, A^(A), and the classical phase function, 4>(A,(a). 
Here, A is the photon wavelength, and a is the exoplanet¬ 
centric angle between the star and the observer. In terms 
of these parameters, the planet/star flux ratio is given by 




( 1 ) 


where Rp is the planet ra dius a nd d is the orbital distance 
(see appendix. Equation ( |A.2[ )). <I>((a, A) is normalized to 
unity at full phase {a = O'"), which defines the geometric 
albedo. The planet’s atmospheric composition, tempera¬ 
ture/pressure structure, cloud/haze profiles, and scatter¬ 
ing properties govern Ag{X) and 4>(A,(a). Hence, obser¬ 
vations of optical planet/star flux ratios, as a function of 


strain the theory of exoplanet atmospheres (e.g., Marley 
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dependent method for measuring orbital parameters that 
complements astrometric and radial-velocity methods. 

In this work, we are interested in the general charac¬ 
ter, as opposed to detailed quantitative predictions, of 
the direct detectability of EGPs. We, therefore, gener- 
ally as sume planets with constant Ag (the exception is in 
^5.1.5, where we vary Ag with orbital distance usi ng two 
simple models) and use the f ormalism presented in Mad- 
husudhan & Burrows (2012) to calculate the associated 
single-scattering albedo [uj) and phase function for an as- 
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Table 1 

Geometric albedos as a function of single-scattering albedo 


(Jj 


A 

9 



Lambert 

Isotropic 

Scalar 

Rayleigh 

Vector 

Rayleigh 

0.005 

0.0033334 

0.00062685 

0.0009394 

0.00094 

0.01 

0.0066668 

0.0012574 

0.0018826 

0.00188 

0.05 

0.033334 

0.0064412 

0.0095705 

0.00963 

0.1 

0.066668 

0.013293 

0.01956 

0.0198 

0.15 

0.1 

0.020604 

0.030015 

0.03057 

0.2 

0.13334 

0.028429 

0.040992 

0.04201 

0.25 

0.16667 

0.036836 

0.052557 

0.0542 

0.3 

0.2 

0.045906 

0.064788 

0.06723 

0.35 

0.23334 

0.055736 

0.077783 

0.08121 

0.4 

0.26667 

0.066447 

0.09166 

0.09628 

0.45 

0.3 

0.078189 

0.10657 

0.11261 

0.5 

0.33334 

0.091154 

0.12269 

0.13041 

0.55 

0.36667 

0.10559 

0.14029 

0.14996 

0.6 

0.40001 

0.12183 

0.15968 

0.17161 

0.65 

0.43334 

0.14033 

0.18131 

0.19585 

0.7 

0.46667 

0.16174 

0.20583 

0.22336 

0.75 

0.50001 

0.18703 

0.2342 

0.25514 

0.8 

0.53334 

0.21777 

0.26799 

0.29281 

0.85 

0.56668 

0.25672 

0.30994 

0.33916 

0.9 

0.60001 

0.30957 

0.36572 

0.39991 

0.95 

0.63334 

0.392 

0.45097 

0.49069 

0.99 

0.66001 

0.53463 

0.59585 

0.64034 

0.999 

0.66601 

0.63582 

0.69783 

0.74337 

1.0 

0.66668 

0.68956 

0.75165 

0.79766 

Note. 

— These calculations were 

carried out with the formal- 


ism described in [Madhusudhan fc Burrows| ( |201^ . 


sumed scattering mechanism. In Table we have tabu¬ 
lated Ag as a function of uj for scalar and vector Rayleigh 
scattering, Lambert reflection, and isotropic scattering. 
Unless stated otherwise, we will hereafter refer to vector 
Rayleigh scattering as Rayleigh scattering. 

The geometric albedo for a Lambert surface has a 
well-known analytic expression: Ag{uj) = (2/3) cj. An¬ 
alytic fits to Ag{uj) for isotropic (iso) and scalar/vector 
(sca/vec) Rayleigh scattering are: 


jiso 

^9 


jsca 

^9 


A 


vec 

9 


0.6896 

0.8190 

0.6901 


(1 +0.4380s) (1 - s) 

(l + 1.2966s)(l +0.7269s)’ 
(H-0.1924s)(l -s) 

(l + 1.5946s)(l + 0.0886w)’ 
(1 +0.4664s) (1-s) 

(l + 1.7095s)(l-0.1345w)’ 


( 2 ) 

( 3 ) 

( 4 ) 


where s = \/l — cu. These fits are accurate to within 0.5% 
for all ay with the largest discrepancies corresponding to 
cj < 0.^ Fits such as these are particularly useful for 
efficient theoretical calculations of the geometric albedo, 
as well as for retrieving a representative atmosphere- 
averaged scattering albedo from an observed geometric 
albedo. 

It is common to assume the Lambert phase function 
when modeling reflected-light observations. However, 
this assumption neglects the albedo- and wavelength- 
dependence of physical phase curves, and generally 
predicts higher values than t he more physically moti 


vated Rayleigh phase function (Madhusudhan & Burrows 


^ [Madhusudhan & Burrows] (|2012| also provide an_analytic fit 
to wJiicJi IS marginally better than Equation Q for small 

cu, but considerably less accurate for oj close to unity. 



Figure 1. Comparison of phase functions for Rayleigh, isotropic, 
and Lambert scattering. The colors correspond to the indicated 
single-scattering albedos (ca). Note that the Lambert phase func¬ 
tion (solid black line) is independent of ca. 

2012|). Figure fl] compares the phase curves for scattering 
due to Rayleigh, isotropic, and Lambert scattering. We 
see that Rayleigh scattering phase curves are systemati¬ 
cally lower than isotropic and Lambert phase curves for 
a < 120°, differing by as much as a factor of ^2 for phase 
angles in the range 6 0° ^ ^ 90° (see also, Madhusud- 

han & Burrows|2012|). For most phase angles, increasing 
cj results in better agreement between these scattering 
processes. Rayleigh phase curves generally increase with 
increasing cj, whereas isotrop ic sca ttering phase curves 
show the opposite trend. In §5.1.1[ we compare predic¬ 
tions of the direct detectability of EGPs with Rayleigh, 
isotopic, and Lambert phase functions. 

2.1.1. Keplerian Elements 

Once ^(A, a) is obtained for a particular planet, the ob¬ 
served phase curve depends on the Keplerian elements of 
its orbifj^ These include the semi-major axis (a), orbital 
inclination (i), eccentricity (e), argument of periastron 
(a;p), time of periastron passage (t^,), and longitude of 
the ascending node (U) (Sndarsky et al. 2QQ5| |Dyndina 


et al. 


scendmg node (bndarsKy et al. 2UUbl [Uyndina 
I2QQ51 [Kane fc Geliiio||2Ull| |lVladhnsudhan ^ 
2012D. For an assumed stellar mass, the period of 


rows 2012p. For an assumed stellar mass, the period 
the planet s orbit (P) can be calculated from a. Note U is 
necessary when an absolute celestial frame is designated 
(e.g., when measuring the full Stokes polarization or cat¬ 
aloguing many different systems in the same solar-system 
frame). However, one can assume Q = 90° or simply ig¬ 
nore this parameter when the system itself is allowed to 
provide a natural orientation (e.g., in exoplanet radial- 
velocity measurements). Throughout this work, we will 
assume D = 90°. 


2.1.2. Albedo Speetra from Homogeneous Atmospheres 

Semi-infinite homogeneous atmospheres are character¬ 
ized by an average single-scattering albedo spectrum: 
a;(A) = Us/{ns + where ns is the scattering opac¬ 
ity, and na is the absorptive opacity. In this section, 

^ See the appendix for a review of the equations that connect the 
Keplerian orbital elements to the observed phase angle and time. 
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Figure 2. Demonstration of the strong dependence of geometric (top row) and scattering (bottom row) albedo spectra upon temperature 
and atmospheric metallicity. The models in the left column fix the temperature at T = 200 K and vary the metallicity, and those in the 
right column fix the metallicity at 10 x solar and vary the temperature. Note the methane features near ~0.62 p,m, ~0.74 /rm, ~0.81 /rm, 
and ~0.89 /rm, ammonia spectral bands near ~0.65 /rm and ~0.79 /rm, and broad water band at ~0.94 /rm. The models are cloud-free, 
homogeneous, serni-infin ite. and use opacities at a pressure of 0.5 atmo spheres. The calcu l ations of cl;(A) use the opacity database from 
[Sharp &Z, Burrows] ( |2007] > and the equilibrium chemical abundances from [Burrows & Sharp] ( 1999|). The geometric albedo spectra assume 
purely KayleigJi scattering atmospheres. 



A [/rm] 


Figure 3. Example variation with wavelength and phase angle of 
the Rayleigh phase function, where we have assumed the single¬ 
scattering albedo spectrum (shown as the black line in the bottom 
left panel of Figure]^ of a semi-infinite homogeneous atmosphere 
with T = 200 K, ^ar metallicity, and a pressure of 0.5 atmo¬ 
spheres. For a Lambert surface, the phase curve does not depend 
upon wavelength, and the horizontal lines corresponding to each 
phase would be fiat and featureless, while of course the correspond¬ 
ing geometric albedos retain a stiff dependence upon wavelength. 


we provide illustrative scattering albedo spectra for cool 


EGPs, which will aid in the interpretation of some of the 
results presented in th is paper. Our calculation s use the 
opacity database from Sharp fc Burroy^ (|2QQ7|) and the 
equilib rium chemical abundances from [Burrows & Sharp] 
( 1 ^. 

Lrom its definition, it is clear that uj{X) is a function 
of atmospheric parameters such as temperature, pres¬ 
sure, and composition. In Figure we demonstrate the 
characteristic strong dependence of scattering albedos 
(bottom row) and the corresponding Rayleigh scatter¬ 
ing geometric albedos (top row) upon temperature and 
metallicity. The models are cloud-free, homogeneous, 
semi-infinite, and use opacities at a pressure of 0.5 at¬ 
mospheres. For the temperatures and pressures relevant 
to cool EGPs, we find ci;(A) is relatively insensitive to 
variations in atmospheric pressure. In general, higher 
temperatures/metallicities result in lower albedos, par¬ 
ticularly for wavelengths longward of ^0.65 /im. 

The spectra in Figure show prominent methane 
features near ^0.62 /im, ^0.74 /im, ^0.81 /im, and 
^0.89 /im, ammonia spectral bands near ^0.65 /im and 
^0.79 /im, and a broad water band at ^0.94 /im. Al¬ 
though there is clearly some degeneracy between metal¬ 
licity and temperature, this figure suggests optical flux 
ratio measurements, particularly between ^0.5 — 0.6 /im 
and/or ^0.82 /im and ^0.94 /im, have great potential as 
temperature/composition diagnostics. 
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(t-tp)/p (t-tp)/p 


Figure 4. The top row shows the time dependence of the planet/star flux ratio for circular orbits with many different inclinations (left) 
and face-on orbits with many different eccentricities (right). The bottom row shows the corresponding planet/star angular separation, 
assuming the system is located 10 pc from Earth. All calculations assume a Jupiter-like planet in an orbit with a = 5 AU and oup = 0°. 


Given a single-scattering albe do spectrum, one can 
use the for malism presented in [Madhusudhan fc Bur-| 
(2012) to calculate the phase function for an as- 


rows 


sumed scattering process as a function of phase angle 
and wavelength. Note that the often-assumed Lambert 
phase function is independent of co’(A), and hence A. In 
Figure we show the phase function for Rayleigh scat¬ 
tering as a function of wavelength and phase angle, where 
we have assumed the scattering albedo spectrum shown 
as the black line in the bottom left panel of Figure 
As expected, the phase curve is equal to one for all A at 
a = 0°, but note the variation in $(A,a) as a increases 
from 0° to 170°. In contrast, a Lambert surface would 
result in featureless, horizontal lines for all a. 


2.2. Planet/Star Angular Separation 

In addition to the planet/star contrast ratio, the 
planet/star angular separation is an essential parameter 
in determining the feasibility of the direct detection of re¬ 
flected light from an exoplanet. For coronagraphic obser¬ 
vations, the figure of merit is given as an IWA and OWA. 
An IWA (OWA) is the minimum (maximum) planet/star 
angular separation that can be resolved by a given instru¬ 
ment. Using the formalism outlined in the appendix, the 
projected angular separation, as a function of phase an¬ 
gle and time, can be obtained via the following equation: 


Sc/) = —— y cos^ {0 ujp) + sin^ {0 + ujp) cos^ i , (5) 

D© ^ 


where D© is the system’s distance from Earth and 0 is the 
planet’s orbital angle as measured from periastron (the 
true anomaly). Inspection of the above equation reveals 
that, as with the planet/star flux ratio, the time depen¬ 
dence of the planet/star angular separation will depend 
sensitively upon the Keplerian orbital elements. 

3. VARIATION OF FLUX RATIOS AND ANGULAR 
SEPARATIONS WITH ORBITAL PARAMETERS 

Using the methods presented in il to calculate 
planet/star flux ratios and angular separations, we now 
present example realizations of the {Fp/F^) — 5(j) param¬ 
eter space for EGPs. The calculations in this section are 
representative of those that will be used in our Monte 
Garlo experiments, which are presented in Q Several 
previous studies ha ve also investigated this py ameter 
space. In particular, Kane & Gelino (2010 2011) consid¬ 
ered the eccentricity an d inclination dependence of exo¬ 
planet phase signatures, |Kane| (|2013|) used the eccentric¬ 
ity distribution of known exopla nets to study t he com¬ 
pleteness of imaging surveys, and Janson (2010) studied 
the effect that prior knowledge of or bital inclination has 
on the efficiency of directly imaging Earth-like planets. 
In this paper, we build on previous work by performing 
a systematic study of the dependence of the direct de¬ 
tectability of EGPs upon a wide range of technological 
and astrophysical parameters (^. 

3.1. Eccentricity and Inclination 
As mentioned in §2.1.1| and ^2.2, observed flux ratios 


and separations are strong functions of the Keplerian or- 
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Figure 5. Planet/star flux ratios (top row) and angular separations (bottom row) for a Jupiter-like planet at a distance = 10 pc from 
Earth on an orbit with e = 0.6 and a = 5 AU. Each column corresponds to orbits with the periastron argument indicated in each top left 
corner, and the colors represent the indicated inclinations. For the orbit orientation, we use the convention cx = 0 = when Up = i = 90°, 
implying full phase corresponds to the time of periastron passage when ujp = 90°. Note that the detectability of planets on face-on orbits 
(purple lines) is independent of ujp. 


bital elements. Figure demonstrates this dependence 
for circular orbits with many different inclinations (left 
column) and face-on orbits with many different eccen¬ 
tricities (right column). The top row shows planet/star 
flux ratios, and the bottom row shows the correspond¬ 
ing angular separations for a system located 10 pc from 
Earth. All calculations in this figure assume a Jupiter- 
like planet in an orbit with a = 5 AU and ujp = 0°. 

Although it has been pointed out previously in the lit¬ 
erature, it is worth reemphasizing that direct-imaging 
observations are well-suited for studying face-on systems 
(the opposite is true for the radial-velocity and transit 
methods, which are biased towards high-inclination sys¬ 
tems). The reason for this is that an exoplanet on a cir¬ 
cular, face-on orbit with FpjF^ > Cmin and > IWA 
will be directly detectable at all times (note we are as¬ 
suming J(/) < OWA ). If, on the other hand, its inclina¬ 
tion is ^90°, both the planet/star flux ratio and angular 
separation will drop below their detection thresholds for 
some fraction of the orbit as the planet passes through 
the full range of possible phase angles. For instance, if 
Gmin = 10“^^ and IWA = 0.3", then a face-on orbit (flat 
purple lines in Figure will be well above the detec¬ 
tion thresholds at all times. However, as we increase the 
inclination in the left column of this figure, we see that 
the flux ratio dips below Cmin for an increasingly larger 
fraction of the orbit centered at new phase {a = 180°). 
Similarly, the minimum planet/star angular separation 
decreases with increasing inclination, falling below the 
IWA for an increasingly larger fraction of the orbit cen¬ 
tered at new phase and full phase. Although the flux 
ratio near full phase is largest for high-inclination or¬ 
bits, the planet will still be undetectable due to the small 
planet/star angular separation. 

In the right column of Figure we show the effect of 
increasing the eccentricity of a mice-on orbit from circu¬ 
lar (e = 0) to highly eccentric (e = 0.8). In this case. 


the phase variations are due to the time dependence of 
the orbital distance. On one hand, a planet on a highly 
eccentric orbit can be an excellent direct-imaging target, 
as it spends a large fraction of its orbit at larger or¬ 
bital distances from its parent star, and for certain sets 
of Keplerian elements, this can provide angular separa¬ 
tions that are highly favorable for direct-imaging obser¬ 
vations. This, however, comes at the cost of lower flux 
ratios throughout most of its orbit and very small angu¬ 
lar separations when it is near periastron passage. 


3.2. Argument of Periastron and Inelination 


As has been pointed out by Kane & Gelino (2011), the 
interaction between the argument of periastron and in¬ 
clination in {Fp/F^) — d(j) parameter space can be quite 
complex. Figure shows various combinations of these 
parameters for a Jupiter-like planet on an eccentric or¬ 
bit with e = 0.6 and a = 5 AU. Each column shows 
planet/star flux ratios and angular separations for orbits 
with inclinations that vary from 0° to 80°, with the argu¬ 
ment of periastron held constant at the value indicated 
in each top left corner. Note we use the convention that 
when a = 0 = t)°^ujp = i = 90°. In other words, perias¬ 
tron passage coincides with full phase when the argument 
of periastron is 90° and with new phase when it is 270°. 

The detectability of planets on face-on orbits is inde¬ 
pendent of cjp, as can be seen by the constant purple 
curve in all panels of Figure As we increase the in¬ 
clination, the flux ratio can either increase or decrease, 
depending upon the periastron argument and phase of 
the orbit. For uOp = 90°, the flux ratio decreases with 
inclination for most of the orbit as the planet transitions 
from being permanently oriented at quarter phase when 
i = 0° to displaying nearly all phase angles when i = 80°. 
At full phase, which corresponds to periastron passage, 
the flux ratio achieves its maximum for this set of pa¬ 
rameters. However, as is typically the case, this maxi- 
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mum flux ratio corresponds to a minimum planet/star 
angular separation that, in the limit of a perfectly edge- 
on system, drops to zero at the time of periastron pas¬ 
sage. Interesting behavior also occurs when uOp = 270°. 
For this periastron argument, periastron passage corre¬ 
sponds to new phase, and, in exactly the opposite sense 
from the uOp = 90° case, the flux ratio increases with in¬ 
clination for most of the orbit and results in a minimum 
at the time of periastron passage for high-inclination or¬ 
bits. For orbits with i > 50°, the minimum angular 
separation actually coincides with the minimum flux ra¬ 
tio, and the maximum angular separation coincides with 
the relatively large fraction of the orbit with a large flux 
ratio. 


4. MONTE CARLO EXPERIMENTS 

We have performed a set of Monte Carlo (MC) experi¬ 
ments, in which we generate random orbits and calculate 
the fraction of time a giant planet on each orbit would be 
directly detectable in the optical. We define this fraction 
to be the observability fraction, /obs, and we calculate its 
value under many different assumptions. An illustrative 
calculation of /obs for a Jupiter-like planet on a random 
orbit around a star located 10 pc from Earth is shown in 
Figure!^ We have assumed a direct-imaging observatory 
with (Xnin = 10“^ and IWA = 0.2", which are the base¬ 
line coronagraph design parameters for WFIRST/AFTA. 
The purple curve represents the planet/star contrast ra¬ 
tio (left axis), the red curve shows the corresponding an¬ 
gular separation (right axis), and the corresponding de¬ 
tection thresholds are indicated by the dashed lines. For 
the planet to be detectable, each curve must be above 
its associated dashed line. This is true only for times 
within the gray hatched region. A Jupiter-like planet on 
this particular orbit would have /obs ~ 20%. Given an 
estimate of the occurrence rate of such planets, the mean 
/obs foi* many systems can be translated into the proba¬ 
bility that a direct-imaging ob serva tion of a random star 
will result in a detection (see ^5.2). 


4.1. Orbital Parameter Distributions 

Each random orbit requires the generation of four Ke- 
plerian elements: {a, e, i, cj^}. To accomplish this, 
we assume planetary orbital orientations are isotropi¬ 
cally distributed with respect to the observer, which 
leads to periastron arguments distributed uniformly in 
[0, 27r] and cosi distributed uniformly in [—1,1]. The 
semi-major axis and eccentricity distributions must be 
determined empirically. Most of our knowled ge of these 
distribution s has come from rad ial-velocity (G umming 
et al. 1 120081 iHoward et al.|l2010|) and transit (|Batalha 
et al.pnl |F ressin et al.||2(J13p surveys that are biased 
towards short-period planets and have only been in oper¬ 
ation long enough to fully describe planets with periods 
up to a few thousand days. 

Nevertheless, we assume both the semi-major axis and 
eccentricity distributions derived from radial-velocity 
planets can be extrapolated to wider separations. For 
the semi-major axis distribution, we adopt a power-law, 
dNjda oc where our fiducial calculations assume 
P = —0.61 (Gumming et al. 2008). This distribution 
is based on radial-velocity planets with orbital periods 
in the range 2 — 2000 days, corresponding to semi-major 



Figure 6. The planet/star flux ratio (purple, left axis) and 
angular separation (red, right axis) of a Jupiter-like planet or¬ 
biting a star 10 pc from Earth on a randomly generated orbit 
with {a, e, i, ujp} = {8.7 AU, 0.14, 71°, 211°}. Assuming 
C'min = 10“^ and IWA = 0.2", which are indicated by the purple 
and red dashed lines, this planet would be directly detectable at 
times within the gray hatched region. The observability fraction 
of this planet is /obs ~ 20%. 


axes in the range 0.03 — 3.0 AU. Results from direct- 
imaging surveys have been used to set model-dependent 
limits on the maximum semi-major axis (umax) to which 
the radial-velocity distribution can be extrapolated. Us¬ 
ing Monte Garlo simulations to examine the null results 
from direct-imag i ng sea rches for EGPs around 118 stars, 
Nielsen & Glose| (|2C)10|) determined Umax ^65 — 234 AU, 
depending on the assumed planet luminosity model. 
Brandt et al. (2014) performed a statistical analysis of 
direct-imaging data from a sample of nearly 250 stars, 
finding Umax ^30 — 100 AU, which again depends on the 
adopted cooling model. With these studies as our moti¬ 
vation, our fiducial calculations assume the semi-major 
axis distribution is given by a power-law with a minimum 
<^min = 0.03 AU and a maximum Umax = 65 AU. 

For the eccentricity distribution, we use the t wo- 
parameter Beta distribution given by [Kipping ( 2013| ): 




( 6 ) 


where T is the Gamma function. Our fiducial calcula¬ 
tions assume a = 1.12 and b = 3.09, which are the long- 
period parameters given by Kipping (2013), who derived 
this distribution by regressing the known cumulative den¬ 
sity function of orbital eccentricities from exoplanets de¬ 
tected from the radial-velocity method. The orbits of 
planets with a < 0.1 AU are expected to be rapidly cir¬ 
cularized due to strong tidal forces. Although we do not 
expect the eccentricities of such planets to influence our 
results, we set e = 0 for any planet with a < 0.1 AU. 

We find the direct detectability of EGPs is highly sen¬ 
sitive to the assumed semi-major axis distribution, but 
relatively insensitive to the eccentricity distribution. Our 
fiducial orbital parameter distributions are summarized 
in Table |2l 


5. RESULTS 
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Figure 7. Distribution of observability fractions (/obs) cis a function of orbital parameters for three illustrative MC experiments, where 
we have assumed Cmin = (left column), Cmin = 5 x 10“^® (middle column), and Cmin = 10“^® (right column). Each experiment 

is composed 5 x 10^ random orbits generated from the distributions given in Table In each panel, dots indicate a single orbit, and the 
color-scale represents fobs: which we calculate for direct observations of a Jupiter-lil^ planet orbiting a star at a distance of = 10 pc 
from Earth, using an instrument with IWA = 0.2", OWA = 1.75", and the minimum contrast indicated at the top of each column. 


Table 2 

Fiducial Orbital Parameter Distributions 


Paramter 

Range 

Distribution 

a 

[0.03 AU, 65 AU] 

dN/da oc 

e 

[0,1) 

Pf}(e;a,b-j^ 

cosi 

[-1, 1] 

uniform 

UJp 

[0, 27r] 

uniform 


Note. — For our fiducial semi-maj or axis distribution, w e 
assume a power-law index /3 = —0.61 ( [Gumming et al^|2008| >. 
For our fiducial eccentricity distribution, we assume a = 1.12 
and b = 3.09, which are the long-period parameters given in 
[Kipping] ( |2013| >. We set e = 0 for any planet with a < 0.1 AU. 

^This_ is the Beta distribution from [Kipping[ ([2013J (Equa¬ 
tion 

5.1. Mean Observability Fraetions 

Here, we present a parametric study of the mean ob¬ 
servability fraction, (/obs), foi* different sets of MC exper¬ 
iments, each composed of 5 x 10^ random orbits. For a 
given population of planets, (/obs) can be interpreted as 
a detection probability, assuming all stars have one such 


planet. In Figurewe show the distribution of fobs as a 
function of orbital parameters for three illustrative MC 
experiments. Each column corresponds to an indepen¬ 
dent experiment, where we have assumed Cmin = 10“^ 
(left column), Cmin = 5 x 10“^^ (middle column), and 
Cmin = 10“^^ (right column). In all three cases, we use 
our fiducial orbital parameter distributions, and calcu¬ 
late fobs foi* direct observations of a Jupiter-like planet 
orbiting a star at a distance of = 10 pc from Earth, 
using an instrument with IWA = 0.2", OWA = 1.75", 
and the minimum contrast specified at the top of each 
column. 

In a given panel, each dot represents a single orbit, 
and the color-scale shows the observability fraction for 
the indicated assumptions about the planet and direct- 
imaging observatory. Eor a fixed minimum contrast, the 
semi-major axis has by far the largest effect on /obs, cre¬ 
ating highly clustered regions of low and high values of 
/obs- The sharp drop in observability interior to 2 AU 
corresponds to an angular separation of 0.2", which is 
the assumed IWA. Note, however, that eccentric orbits 
with a < 2 AU are able to achieve nonzero observability 
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fractions, since a planet on such an orbit spends a large 
fraction of its time near orbital distances of a(l + e), 
making it possible for the observed angular separations 
to fall outside the IWA, provided the full set of Keple- 
rian elements allows for such geometries. These results 
make it clear that /obs is very sensitive to the assumed 
semi-major axis distribution. Steeper (shallower) power- 
law indices will result in narrower (broader) regions of 
high /obs values, and the assumed values of Umax and 
amin determine the region within which the 5 x 10^ plan¬ 
ets are distributed. In addition, decreasing the minimum 
contrast increases the range of semi-major axes that are 
detectable and leads to an increasing number of orbits 
for which /obs^1-0 (be., orbits that are observable at all 
times). For Cmin = 10“^, low- to modest-eccentricity or¬ 
bits with /obs >1% are confined to 2 AU ^ ^ 10 AU. 

As we decrease Cmin from left to right in this figure, 
we see orbits with increasingly larger semi-major axes 
achieve nonzero /obs, reaching a ^ 15 AU and a ^ 35 AU 
for Cmin = 5 X 10“^^ and Cmin = 10“^^, respectively. 

As mentioned previously, low-inclination orbits are 
beneficial to direct-imaging observations. In the top left 
panel, it is particularly clear that orbits with /obs ^1-0 
are associated with low inclinations. However, planets 
with i > 50° can have /obs ^15% out to a ^10 — 30 AU, 
depending on Cmin, since phase angles near full phase 
become observationally accessible. There is also an inter¬ 
esting feature near ^270°, which corresponds to orbits 
with periastron passage occurring at new phase. Near 
this periastron argument, the region of large /obs extends 
to semi-major axes as large as ^20 AU for Cmin = 10“^^. 
For a planet with ^ 270°, the planet/star flux ra¬ 
tio rises with increasing inclination throughout most of 
its orbit, and the largest angular separations coincide 
with the highest flux ratios when i > 50° (see Figure]^. 
These factors collectively increase the semi-major axes 
for which a planet on such an orbit can have a large 
observability fraction. 

For the remainder of this section, we present the mean 
observability fraction, (/obs), for many MC experiments, 
which are run under different sets of assumptions. In 
each case, our goal is to study how the observability 
fraction depends upon input parameters such as albedo, 
planetary radius, distance from Earth, minimum achiev¬ 
able contrast, and the semi-major axis distribution. For 
each experiment, we vary one parameter with respect to 
a fiducial model, which is comprised of the following as¬ 
sumptions: 

• the planetary atmospheres are Rayleigh scattering 
with Ag = 0.5 (this corresponds to cj ~ 0.95 for 
Rayleigh scattering), 

• the planets have the radius of Jupiter and orbit 
stars located 10 pc from Earth, 

• the direct-imaging observatory has an OWA = 
1.75" and Cmin = 10“^ (similar to what is antici¬ 
pated for WFIRST/AFTA), 

• the Keplerian orbital elements are distributed ac¬ 
cording to the distributions given in Table 

Of course, the above assumptions do not fully represent 
realistic physical/statistical properties of planets or any 


future direct-imaging mission. Nevertheless, by varying 
each parameter of this very simple reference model, it 
is possible to gain insight into the general problem of 
directly detecting giant exoplanets in the optical, both 
in the case of a blind search and observations of stars 
known to host giant planets. 

5 . 1 . 1 . Dependence on Scattering Properties 

We start by varying the first of our fiducial model 
assumptions. Figure shows the dependence of 
(/obs) upon the geometric albedo (left panel) and the 
single-scattering albedo (right panel), assuming Lam¬ 
bert (black), isotropic (red), and Rayleigh (blue) scat¬ 
tering atmospheres. Each curve is composed of many 
distinct MC experiments, and the three sets of curves 
in each panel correspond to the indicated IWAs. The 
WFIRST/AFTA coronagraph is anticipated to have 
IWA ^ 0.2". Note this value will depend upon the an¬ 
gular separation and wavelength of the observation. For 
each scattering mechanism, there is a one-to-one corre¬ 
spondence between uj and Ag (see TablefB. Additionally, 
the phase curves for isotropic and Raylei^ scattering are 
direct functions of cj, and the Lambert phase curve is in¬ 
dependent of UJ (see Figure!^. Thus, given Ag or cj, 
it is straightforward to calculate the product Ag • ^((a), 
which is proportional to the planet/star flux ratio, for an 
assumed scattering or effective scattering mechanism. 

In the left panel of Figure we show the dependence 
of (/obs) upon Ag. For each MC experiment, we assume 
Ag and interpolate Tabled to infer the corresponding cj, 
which is needed to calculate the phase curves for Rayleigh 
and isotropic scattering. For a given A^, the scatter¬ 
ing phase functions are entirely responsible for the dif¬ 
ferences seen between the assumed scattering processes. 
Note that, since the Lambert phase function is indepen¬ 
dent of cj, the phase function is the same for each black 
curve. As can be seen in Figure Rayleigh scattering 
phase curves are systematically lower than isotropic and 
Lambert phase curves for a < 120°. It is for this rea¬ 
son that (/obs) is always the lowest for Rayleigh scatter¬ 
ing. For most A^, Rayleigh scattering predicts observ¬ 
ability fractions that are ^20% below the predictions of 
isotropic and Lambert scatting, regardless of the IWA. 
At the anticipated IWA for WFIRST/AFTA, (/obs) de¬ 
creases sublinearly with Ag for each scattering process, 
with maximum and minimum values for Rayleigh scat¬ 
tering of ^11% for Ag = 0.65 and ^1% for Ag < 0.1. For 
context. Figure ^suggests that cool EGPs may have large 
geometric albedos (A^^ 0.6 —0.8) in the optical, particu¬ 
larly for wavelengths ^0.6 /im and near the broad water 
band at ^0.94 /im. As expected, these results suggest 
planetary geometric albedos can have a dramatic effect 
on the results of an optical direct-imaging survey. 

In the right panel of Figure we show the dependence 
of (/obs) upon UJ. In contrast to the left panel, here we 
assume uj for each MG experiment and use Table to 
infer the corresponding Ag. In this case, both Ag and 
4>((a) contribute to the differences between the assumed 
scattering processes. We see that (/obs) is comparable for 
Rayleigh and isotropic scattering for all uj. This arises 
from the fact that, for a given cj, Ag for Rayleigh scatter¬ 
ing is systematically higher than that for isotropic scat¬ 
tering, whereas 4>(a) for Rayleigh scattering is system¬ 
atically lower^ leading to roughly consistent observabil- 
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Figure 8. The dependence of the mean observability fraction (fobs) upon the geometric albedo (Ag, left) and single-scattering albedo 
(ca, right), assumir^ Lambert (black), isotropic (red), and Rayleigh (blue) scattering atmospheres. In the left panel, we assume Ag and 
interpolate Table [Jto infer the corresponding u for each scattering process. In contrast, in the right panel, we assume cj and infer Ag for 
each scattering process. For each calculation, we assume Jwiter-size planets that orbit stars located 10 pc from Earth, with random orbital 
parameters drawn from our fiducial distributions (Table 1^. The direct-imaging observatory is assumed to have an OWA = 1.75'' and 
C'min = 10“^. The three sets of curves in both panels correspond to the indicated IWAs. The WFIRST/AFTA coronagraph is anticipated 
to have IWA = 0.2". 



Figure 9. Variation of the mean observability fraction (fobs) 
with planet radius. With the exception of p lanet radius, all as¬ 
sumptions are as listed at the end of §5.I| The three curves 
show (fobs) for observations made with a coronagraph with the 
specified IWAs. We note that gas giants are expected to have 
0.8 Rjup ^ ^ 1-2 Rjup: but we show a wider range of radii for 

illustrative purposes. 

ity fractions. For conservative scattering {uj = 1.0), the 
geometric albedos for Lambert, isotropic, and Rayleigh 
scattering differ by less than ^13%, and {fobs) is essen¬ 
tially independent of the assumed scattering mechanism. 
However, for non-conservative scattering {uj < 1.0), the 
geometric albedos for isotropic and Rayleigh scattering 
can be less than those for Lambert scattering by as much 
as a factor of ^4. As a consequence, Lambert scattering 
generally predicts significantly higher (/obs) than both 
isotropic and Rayleigh scattering. 

As a specific example, consider uj = 0.6. For this 
scattering albedo, the Lambert phase function corre¬ 


sponds to a geometric albedo that is a factor of ^3 
larger than isotropic scattering and a factor of ^2 larger 
than Rayleigh scattering. At the WFIRST/AFTA IWA 
of 0.2", this results in (/obs) ^ 3% for Rayleigh and 
isotropic scattering and (/obs) ^ 10% for Lambert scat¬ 
tering. Thus, when the geometric albedo is calculated 
self-consistently with the single-scattering albedo, the as¬ 
sumed scattering phase function can greatly influence the 
predicted yield from an optical direct-imaging survey. In 
particular, assuming the Lambert phase function gen¬ 
erally results in an overestimate of the detectability of 
giant planets with respect to the more physically moti¬ 
vated Rayleigh phase function. 

5 . 1 . 2 . Dependence on Planet Radius 

The planet/star flux ratio is proportional to the square 
of the planet radius. Therefore, we expect (/obs) fo 
be sensitive to small variations in planet size, which we 
have heretofore assumed to be that of Jupiter. Figure 
shows (/obs) a.s a function of planet radius, assuming ob¬ 
servations made with a coronagraph that has IWA = 
O.I", 0.2", and 0.3". All assumptions other than planet 
radius are as listed at the end of §5.I[ We note that gas 
giants are expected to have 0.8 Rjup ^ Rp ^ 1-2 Rjup, 
but we show a wider range of radii for illustrative pur¬ 
poses. Within the range of expected radii for gas giants, 
an increase in radius by a factor of ^1.5 leads to an in¬ 
crease in (/obs) by a factor of ^1.5 — 2.5, depending on 
the IWA. We find that smaller IWAs are less sensitive 
to variations in planet radius. This effect likely arises 
from the d~‘^ dependence of the planet/star flux ratio, 
which—because of the increasing number of planets at 
smaller orbital distances that are accessible to smaller 
IWAs—reduces the weight carried by the dependence 
of the flux ratio. 

5 . 1 . 3 . Dependence on Distance and Minimum Contrast 
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Figure 10. The mean observability fraction (fobs) cis a function 
of distance from Earth (D^) and minimum achievable contrast 
(Dmin), assuming the IWA = 0.1", 0.2", and 0.3". Dashed lines 
assume Cmin = 10“^, and solid lines assume Cm i n = 10“^° . All 
other assumptions are as listed at the end of §5T| 


The fiducial model listed at the end of ^5.1 assumes 


direct observations of planet/star systems located a dis¬ 
tance = 10 pc from Earth with a minimum achiev¬ 
able contrast Cmin = 10“^. It is obvious that the di¬ 
rect detectability of EGPs will vary strongly with both 
of these parameters. In Eigure we investigate this 
variation in the context of the mean observability frac¬ 
tion. We show (/obs) foi* planet/star systems located 
0.01 — 50 pc from Earth, assuming minimum achiev¬ 
able contrasts of 10“^ (dashed lines) and 10“^^ (solid 
lines). The calculations assume a coronagraphic instru¬ 
ment with an IWA = 0.1", 0.2", and 0.3". All other 
parameters are as listed at the end of §5.1[ 
WEIRST/AETA is anticipated to have IWA ^0.2" and 
Gmin ^ 10“^. In this case. Jupiter-like planets located 
10 pc from Earth spend on average ^10% of their time 
in observable configurations, whereas the same planets 
located 30 pc from Earth spend on average ^0.1% of 
their time in observable configurations. If instead a min¬ 
imum contrast of 10“^^ can be achieved, then these frac¬ 
tions become ^30% and ^16%, respectively. At small 
distances from Earth, we see (/obs) increases with dis¬ 
tance. This trend is the result of planets at sufficiently 
wide orbital distances to spend some fraction of their 
orbit outside the OWA, which we have assumed to be 
1.75". Assuming a uniform distribution of stars and 
IWA = 0.2", Jupiter-like planets orbiting stars within 10, 
30, and 50 pc from Earth have volume-averaged observ¬ 
ability fractions of ^12%, 3%, and 0.5% for Cmin^lO”^ 
and ^28%, 20%, and 13% for 


5.1.4. Dependence on the Semi-Major Axis Distribution 

Our fiducial model assumes the semi-major axis dis¬ 
tribution derived from radial-velocity planets can be ex¬ 
trapolated to wider separations (a > 3 AU). However, 
the validity of this assumption is highly uncertain, and, 
as evidenced by Eigure the semi-major axis distri¬ 
bution has a large influence on the distribution of ob¬ 
servability fractions. In this section, we explore the de¬ 
pendence of (/obs) upon the semi-major axis distribution 


by varying each of its components: the minimum semi¬ 
major axis (amin), the maximum semi-major axis (amax), 
and the power-law index of the distribution 

Our fiducial model assumes amin = 0.03 AU (note 
we set e = 0 for any planet with a < 0.1 AU), which 
has been in ferred for radial-velocity planets (|Cumming| 
et al. |2008). In Eigure we show (/obs) as a func- 


tion of UmiL assuming the indicated IWAs. At a given 
IWA, the mean observability increases with increasing 
amin until Umin/J^© ^ IWA, whcrc = 10 pc in our 
fiducial model. Since all of the MG experiments use the 
same number of planets (5 x 10^), increasing the min¬ 
imum semi-major axis forces more planets to exist at 
larger orbital separations, leading to larger (/obs) when 
Umin/^© < IWA and smaller (/obs) when amin/D^ > 
IWA. At the anticipated IWA of WEIRST/AETA, the 
mean observability varies by ^1% for Umin iu the range 
0.03 - 2.0 AU. 



of planets, smaller Umax values yield larger {fobs) val¬ 
ues, which is sensible given that more planets will ex¬ 
ist at smaller, but still detectable, orbital distances. As 
Umax increases, (/obs) slowly decreases as a result of the 
increasing number of planets at very large orbital dis¬ 
tances, which are too faint to be detected and, in the 
case of nearly face-on orbits, fall outside the OWA. Over 
the range Umax ^ 30 — 100 AU, (/obs) decreases by a 
factor of ^1.6 for each IWA. 

The power-law index of the semi-major axis distribu¬ 
tion of radial-velocity planets with periods in the range 
2 — 2000 days has been constrained to be /3 = —0.61 ±0.15 
( Gumming et al.||2008 ). Our fiducial model assumes this 
index can be extrapolated to planets at much larger or¬ 
bital separations. In Eigure we test the sensitivity 
of (/obs) to variations in the power-law index. As we 
steepen P towards more negative values, an increasing 
number of planets begin to pile up at small orbital sep- 



Figure 11. Dependence of the mean observability fraction (fobs) 
upon the minimum semi-major axis (cimin) of fho semi-major axis 
dis tribu tion. All other model parameters are as listed at the end 

of O 
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Figure 12. Dependence of the mean observability fraction (fobs) 
upon the maximum semi-major axis (a,max) of the semi-major axis 
dis tribu tion. All other model parameters are as listed at the end 
of ^ 



p 


Figure 13. Dependence of the mean observability fraction (fobs) 
upon the power-law index of the semi-major axis distrib ution (/3). 
All other model assumptions are as listed at the end of O 

arations, which leads to larger {fobs) values. At the 



/3, the difference between the minimum and maximum 
{fobs) values is only ^3%. 


5.1.5. Geometric Albedo as a Function of Distance 

Thus far, we have assumed the geometric albedo is con¬ 
stant throughout the orbit of each planet. However, the 
atmospheric temperature of a planet varies strongly with 
orbital distance, and such variations determine whether 
water or ammonia clouds will form. We, therefore, ex¬ 
pect the albe do spectrum to vary s ignificantly with or¬ 
bital distance (Sudarsky et al.|2005). We also expect the 
albedo spectrum to depend sensitively upon the planet’s 
mass, age, and parent star. In this section, we compare 


our fiducial assumption of = 0.5 = constant with two 
simple models that account for the distance dependence 
of the albedo spectrum. 

The first mo del is an interpolation o f the models shown 
in Figure 9 of Sudarsky et al.| (2005), which depicts the 
geometric albedo spectra of a Jupiter-mass, 5-Gyr planet 
orbiting a G2 V star at distances in the range 0.2 — 15 AU. 
In Figure [T^ we show the predictions of this model for 
the distance dependence of the geometric albedo for var¬ 
ious wavelengths, including wavelengths associated with 
methane (0.74 /im and 0.89 yum), ammonia (0.65 /im), 
and water (0.94 /im). We see the geometric albedo, 
and hence the planet/star flux ratio, is a strong, non¬ 
monotonic function of both distance and wavelength. 

Using this model for we carried out MG experi- 
ments for observations at the wavelengths shown in Fig¬ 
keeping all other parameters as listed at the end 
For orbital distances d > 15 AU, we set Ag 


of ^5 
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equal to its value at 15 AU; this does not significantly 
influence our results, since planets at these separations 
contribute very little to {fobs) (see Figure Q. In Fig- 
ure we show (/obs) as a function of IWA assuming 
the ^^udarsky et al.” albedo model (solid lines) and 
our fiducial assumption of Ag = 0.5 = constant (solid 
black line). As can be seen in the first column of Fig¬ 
ure most planets with nonzero /obs values have semi¬ 
major axes in the range 2 — 10 AU, so it is the variabil¬ 
ity of Ag in this region that has the strongest effect on 
(/obs)- In the wavelength range we have considered, this 
model suggests the optimal observability fractions occur 
for A ^ 0.5 — 0.6 /im, which Figure indicates is the 
effective bandpass of our fiducial assumption Ag = 0.5. 
At the IWA of WFIRST/AFTA, direct observations in 
the wavelength range A ^ 0.5 — 0.6 /im of a Jupiter- 
like planet located 10 pc from Earth will, on average, be 
possible for ^10% of its orbital period. We note that 
the observed reflected fluxes shortward of ^0.6 /im from 
Jupiter and Saturn are less than what this model would 
predict, which is made man ifest by the reddish hues of 


these solar-system planets (Karkoschka 1994 Burrows 


2014). 


Tor the second model of the distance dependence of the 
albedo spectrum, we calculate the equilibrium tempera¬ 
ture for each planet as a function of orbital distance in 
a self-consistent fashion with the albedo spectra. This is 
accomplished as follows. For each planet in an MG exper¬ 
iment, we assume the planetary atmosphere is homoge¬ 
neous and Rayleigh scattering. In addition, solar metal- 
licity and a pressure of 0.5 bar are assumed (albedo spec¬ 
tra for such a planet with T = 200 K are shown as black 
lines in the left column of Figure j^. At a given phase of 
the planet’s orbit, we assume an initial Bond albedo and 
calculate the equilibrium temperature. Assuming this to 
be the temperature of the atmosphere, we generate (scat¬ 
tering, geometric, and sp herica l) albedo spectra with 
the methods described in §2.1.2[ We then integrate the 
spherical albedo spectrum to obtain the Bond albedo and 
calculate a new equilibrium temperature (Teq). We iter¬ 
ate this process until we have achieved self-consistency 
between the equilibrium temperature and albedo spec¬ 
tra. This procedure is performed at every phase of the 
planet’s orbit, resulting in a geometric albedo spectrum 
and phase function that varies with orbital distance. The 
results of MG experiments assuming this “Teq model” 



































The Direct Detectability of Giant Exoplanets in the Optical 


13 



ci[AU] 


Figure 14. The geometric albedo as a function of distance and 
wavelength for a Jupiter-mass, 5-Gyr planet orbiting a G2 V star. 
T hese results are an inte rpolation of the models shown in Figure 9 
of |Sudarsky et ^ ( |2005] >. 



IWA [arcsec] 


Figure 15. Gomparison of (fobs) assuming three models for the 
albedo spectrum. The “Sudarsky et al.” (solid lines) and “Tgq” 
(dashed lines) models account for the distance dependence of the 
albedo spectrum, whereas the Ag = 0.5 model (solid black line), 
which is our fiducial model, holds it constant throughout each 
planet’s orbit. The solid cyan and purple lines are overlapping, 
which indicates that orbital distances d < 4 AU are d omi nating 
the mean observability at these wavelengths (see Figure [T^. 

are shown as dashed-lines in Figure We see that 
this model is roughly consistent with the Sudarsky et ah 
model for A < 0.7 /im, but its predictions are very differ¬ 
ent at longer wavelengths. This model also predicts that 
the wavelength range A ^ 0.5 — 0.6 /im offers the optimal 
observability fractions. 

5.2. Detection Probabilities 

Given assumptions about the orbital parameter distri¬ 
butions and occurrence rates (fp) of a population of plan¬ 
ets, the mean observability fractions presented in this 
work can be translated into the probability that a blind 
search for such planets will yield a detection. Statistical 


studies of planets discovered by the radial-velocity and 
microlensing methods suggest that the frequency of giant 
planets with a < 20 AU falls in the range fr, ^ 10 — 20% 
(Gumming et al. 2008 Gassan et al. 2012). If we as- 
sume fp = 157o, then calculations based on our fidu- 
cial model assumptions, which are listed at the end of 
^5.1, suggest that the probability of detecting a Jupiter- 
like planet in a blind search around a star located 10 pc 
from Earth is p = fp • {fobs) ^ 1-5%, where we have as¬ 
sumed IWA = 0.2". If a minimum contrast of 10“^^ is 
assumed, then this probability increases to p ^ 4.2%. If 
the star is located 20 pc from Earth, we find p ^ 0.4% 
for Cmin = 10“^ and p = 3.2% for Cmin = 10“^^. The 
scaling of p with parameters such as the IWA, planet 
radius, distance from Earth, and atmospheric scattering 
properties can be inferred from the figures presented in 
this paper. In all but the most optimistic configurations, 
the probability for detection in a blind search carried out 
with WFIRST/AFTA’s baseline coronagraphic capabili¬ 
ties is low (< 5%). 

By the time direct detection in the optical becomes 
feasible, long-term radial-velocity campaigns will have 
likely discovered many wide-separation EGPs, for which 
a subset of Keplerian elements will be known. Further, 
it has been estimated that Gaia will discover and deter¬ 
mine robust orbits for ^1 500 planets with period s < 6 yr 
within 50 pc from Earth (Perryman et al.|2014), includ¬ 
ing ^100 gian t planets around know n IVl dwarf host stars 
within 30 pc (Sozzetti et al. 2014). Direct-imaging sur¬ 
veys that leverage this intorrnation can improve their de¬ 
tection probabilities with respect to blind searches, and 
the tools we have developed for this work can be used 
to determine both the most promising systems to target 
and when to observe them. It is important to point out 
that, as can be inferred from Figuresthe window of 
time that an EGP is directly detectaWe can be very brief 
and sparsely distributed throughout its orbit, suggesting 
that previous orbital parameter constraints must be very 
precise if they are to be used to optimize direct-imaging 
surveys. For a detailed study of the impact of uncer¬ 
tain radial-velocity orbital parameters on direct i maging 
and sp ectroscopy of radial-velocity exoplanets, see|Brown| 

(I2M^. 


5.3. Thermal Emission from Young Planets 

Throughout this work, we have focused on cool, rela¬ 
tively old giant planets, for which the flux from thermal 
emission is negligible compared to the flux from reflected 
starlight. The thermal emission from a giant planet is 
age, mass, and metallicity dependent. Furthermore, the 
particular formation mechanism of giant plane ts signifi¬ 
cantly influences the cooling of the planet (e.g.. Burrows 


W _ _ _ 

et al.|1997l|Marley et al. [20(171 |Fortney et al.|20U8[|bpieg'^ 


6 Burrows 2012|). Regardless of the assumed cooling 

model, it is reasonable to expect the thermal emission 
contribution to the optical planet/star flux ratio will in¬ 
crease with planet mass and decrease with planet age. 

In Figure we show the planet/star flux ratio 
(top panel) ana associated Rayleigh scattering geomet¬ 
ric albedo spectrum (bottom panel) for a representative 

7 Jupiter-mass planet orbiting a Sun-like star at an or¬ 
bital distance of 5 AU, where we have included the effect 
of blackbody thermal emission from the planet. We use 
the evolutionary calculations of Burrows et al. (1997) to 
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Figure 16. The planet/star flux ratio (top) and geometric albedo 
spectrum (bottom), as a function of wavelength and age, for a 7 
Jupiter-mass planet orbiting a Sun-like star at an orbital distance 
of 5 AU, where we have included the effect of blackbody thermal 
emission from the planet. The black dashed lines in the top panel 
show the component of the contrast ratio that is associated with 
thermal e mission from the plan et. We use the evolutionary calcu¬ 
lations of|Burrows et al.|(|1997|) to model the evolution of the age, 
radius, arid temperature ot tire planet. The atmospheric pressure 
is 0.5 bar, and solar metallically is assumed. We have assumed 
both the star and planet are uniformly luminous blackbodies. The 
albedo spectrum was ca lculated with the opacity database from 
I Shar p fc Burrows (|2Q07|> and t he equilibrium chemical abundances 
trom |l^urrows ^ToJiarpI (|1999|). The atmosphere is assumed to be 
Rayleigii scattering. i'Jie planet age and effective temperature (Tp) 
color-scales apply to both panels. 


model the evolution of the age, radius, and temperature 
of the planet, and the albedo spectrum was calculated 
with the opacity database from Sharp fc Burrows (|2QQ7 ) 
and equilibri um chemical abundances from [Burrows &:| 


Sharp (1999). We assume solar metallically and an at- 
mospheric pressure of 0.5 bar. Both the planet and star 
are assumed to be uniformly luminous blackbodies. That 
is, we assume the thermal emission component of the con¬ 
trast ratio is given by RpBx{Tp)/R‘IB x{T^)^ where Tp is 


the effective temperature of the planet and Bx{T) is the 
Planck function. 

In this simple model, the thermal emission from the 
planet noticeably contributes to the optical planet/star 
flux ratio only for ages < 80 Myr. Nonetheless, we see it 
briefly dominates the flux ratio at wavelengths ^0.5 /im 
when the planet is a mere few tens of millions of years 


old, resulting in a planet that transitions from red to 
blue with increasing age and decreasing effective tem¬ 
perature. The age range over which thermal emission is 
important will be larger (smaller) for giant planets more 
(less) massive than 7 Jupiter masses. The decrease in flux 
at short wavelengths (^0.4 —0.5 jam) for ages > 100 Myr 
is due primarily to the contraction of the planet’s radius, 
which decreases from ^1.4 Rjup to ^1 Rjup throughout 
its evolution. In contrast, the geometric albedo spectrum 
is approximately constant over this same wavelength and 
age interval. 

6. SUMMARY AND CONCLUSIONS 

With the planning of future exoplanet-focused missions 
such as WFIRST/AFTA well underway, it seems likely 
that the direct detection of EGPs in the optical will be¬ 
come feasible in the coming decade. Therefore, studies 
of the observational signatures and direct detectability 
of reflected light from EGPs are particularly timely and 
will aid in the design and optimization of these missions. 
In this work, we have presented an exploration of the di¬ 
rect detectability of EGPs at visible wavelengths, quan¬ 
tified by the fraction of a planet’s orbit for which it is 
in an observable configuration—its observability fraction 
/obs- Using a suite of Monte Garlo experiments, we in¬ 
vestigated the dependence of /obs upon various techno¬ 
logical and astrophysical parameters such as the inner 
working angle and minimum achievable contrast of the 
direct-imaging observatory; the planet’s scattering phase 
function, geometric albedo, single-scattering albedo, ra¬ 
dius, and distance from Earth; and the semi-major axis 
distribution of EGPs. Unless stated otherwise, we now 
summarize our results assuming observations of Jupiter- 
like planet^ using an optical coronagraphic instrument 
with IWA = 0.2" and Cmin = 10“^, which are the base¬ 
line design parameters for WEIRST/AETA. 

Low-inclination orbits are highly favorable for direct- 
imaging observations (Eigures fflandl^. Howeverjplan- 
ets with i > 50° and uOp ^ 2T& (Eigures and are 
generally observable out to larger semi-major axes than 
planets in other configurations. At a distance of 10 pc 
from Earth, only (low- to modest-eccentricity) EGPs 
with semi-major axes in the range a ^ 2—10 AU will have 
fobs > !%• If Gmin ^ 10“^^ Can be achieved, this range 
extends out to a ^ 35 AU, with a much larger fraction of 
orbital parameter space yielding /obs ^ 1 (EigureJ^. In¬ 
spection of the bottom row of Eigurej^ suggests that the 
mean direct detectability of EGPs is rdatively insensitive 
to their eccentricity distribution. 

The mean observability fraction, (/obs), of EGPs is a 
strong function of their atmospheric scattering proper¬ 
ties. Eor a given geometric albedo, Rayleigh scatter¬ 
ing atmospheres predict (/obs) values that are gener¬ 
ally r\j 20% below the predictions of isotropic and Lam¬ 
bert scattering atmospheres (left panel, Eigure 1^. If 
instead one assumes a single-scattering albedo and self- 
consistently calculates the phase function for a given 
scattering mechanism, the often-assumed Lambertian 
phase function can overestimate the detectability of 
EGPs by as much as a factor of ^4 with respect to the 

^ We remind the reader that we deflne a “Jupiter-like” planet to 
be an EGP with the radius of Jupiter, Ag = 0.5, and a Rayleigh 
scattering atmosphere. 
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more physical^ motivated Rayleigh phase function (right 
panel, Figure]^. 

Within the range of expected radii for gas giants 
(O.SRjup ^ Rp ^ l-2Rjup), an increase in radius by 
a factor of ^1.5 leads to an increase in {fobs) by a factor 
of ^1.8 (Figure]^. We find that smaller IWAs are less 
sensitive to variations in planet radius. This effect likely 
arises from the d~‘^ dependence of the planet/star flux 
ratio, which—because of the increasing number of plan¬ 
ets at smaller separations that are accessible to smaller 
IWAs — reduces the weight carried by the dependence 


of the flux ratio. 

An EGP located 10 pc from Earth will spend on av¬ 
erage ^10% of its time in observable configurations, 
whereas the same planet located 30 pc from Earth will 
spend on average ^0.1% of its time in observable con¬ 
figurations. If a minimum contrast of 10“^^ can be 
achieved, then these fractions become ^30% and ^16%, 
respectively (Figure 10). Assuming a uniform distribu¬ 
tion of stars, EGPs orbiting stars within 10, 30, and 
50 pc from Earth have volume-averaged observability 
fractions of ^12%, 3%, and 0.5% for Cmin ^ 10“^ and 
-28%, 20%, and 13% for Cmin - 10“^^ Thus, the sci¬ 
entific return from a future optical direct-imaging survey 
will depend critically upon its minimum achievable con¬ 
trast, as well as the number of wide-separation EGPs 
that exist w ithin a few tens of parsecs from Earth. 

In ^ 5.1.4 we explored the dependence of (fobs) upon the 
semi-major axis distribution of EGPs by varying each of 
its components: the minimum semi-major axis (amin), 
the maximum semi-major axis (amax), and the power- 
law index of the distribution {P). Given constraints from 
radial-velocity planets and infrared direct-imaging sur¬ 
veys, the uncertainty in (fobs) is dominated by the un¬ 
certainty in amax- For Umax iu the range 30 — 100 AU, 
which is the model-dependent limit set by [Brandt et al.| 
(2014|), {fobs) varies by a factor of —1.6 (Figure |l2|). 

4'he albedo spectrum of a giant planet is expected 
to be a non-monotonic fu nction of both distanc e and 
wavelength (Figure and |Sudarsky et al.||2005|). Us¬ 
ing two simple modSsior the distance- and wavelength- 
dependence of the geometric albedo, we find that bluer 
wavelengths result in larger {fobs) values, which is to 
be expected given that we are interested in reflected- 
light observations. We caution that, as revealed by their 
reddish hues, the observed reflected fluxes shortward of 
—0.6 /im from Jupiter and Saturn are less than what our 
models would predi ct. 

As described in §5.2[ the mean observability fraction 
for a population of planets is proportional to the proba¬ 
bility (p) that a blind search for such planets will yield 
a detection: p = /^ • {fobs)^ where fp is the planet oc¬ 
currence rate. The scaling of p with parameters such as 
IWA, Cmin, planet radius, distance from Earth, and at¬ 
mospheric scattering properties can be inferred from the 
figures presented in this paper. In all but the most op¬ 
timistic configurations, the probability for detection in a 
blind search carried out with WFIRST/AFTA’s baseline 
coronagraphic capabilities is low (< 5%). 

Finally, we note that thermal emission from young 
EGPs can provide a non-negligible contribution to op¬ 
tical planet/st ar flux ratios. Assumi ng the evolutionary 
calculations of [Burrows et ah' (1997), the thermal emis¬ 


sion from a 7 Jupiter-mass planet orbiting a Sun-like star 
briefly dominates the flux ratio at wavelengths ^0.5 pm 
when the planet is a few tens of millions of years old, be¬ 
coming completely negligible after —80 Myr. The result 
is a planet that transitions from red to blue with increas¬ 
ing age and decreasing effective temperature (Figure 16). 
The age range over which the thermal emission contribu¬ 
tion to the optical flux ratio is significant will be larger 
(smaller) for giant planets more (less) massive than the 
representative 7 Jupiter-mass planet we have considered 
here. 
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APPENDIX 


The orbital phase angle (a) is related to the true anomaly (0) and the Keplerian elements 

cos a = sin(^ + ujp) sin i sin Cl — cos Cl cos{0 + cjp), 

= sin(^ + Up) sini, 


(AT) 


where i is the orbital inclinati on, u ^^ is the argument of periastron, and Cl is the longitude of the ascending node. The 
second equality in Equation (A.l) comes from the assumption that Cl = 90°. For the orbit orientation, we use the 
convention a = 0 = 0° when Up = i = 90°. The range of observable phase angles depends only on the inclination 
of the orbit: 90° — i < a < 90° + i. For eccentric orbits, the light curve will vary even when i = 0°. This is due to 
changes in the orbital distance (d) along the orbit, which is given by 


d = 


a(l — e^) 
1 + e cos 0 


(A.2) 


The Keplerian elements are connected to time (t) through the mean anomaly (M = 27r(t — tp)/P, where tp is the time 
of periastron passage), which is related to the eccentric anomaly (E) by Kepler’s equation: 


M = E — esinE. 

E can be expressed in terms of the true anomaly as 


sinP 


We can then relate the planet’s true anomaly to time: 


sin^Vl — ' 


1 + e cos 0 


(A.3) 

(A.4) 


t tr) 


1 

27r 


esm 


in^Vl — ' 


1 + e cos 0 


+ 2 tan 


/(1-e) 0 

) -f tan- 

(1 + e) 2 , 


(A.5) 


Thus, for an arbitrary orbit orientation. Equations (|A.1|) and (|A.5|) combine to yield the exact orbital phase at any 
time. 


^ Note D is necessary when an absolute celestial frame is des¬ 
ignated. However, one can assume D = 90° or simply ignore this 


parameter when the system itself is allowed to provide a natural 
orientation. 

















